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SUMMARY 

A study has been made of the ejffects of two nonlinear stability 
derivatives j caused by the nonlinear varixitions of pitching- 
moment and lift coefficients with angle of attack j on the longitudi- 
nal motions of an aircraft. Theoretical methods involving the 
Laplace transformation and the procedures of nonlinear 
mechanics have been presented along with electrical-analog re- 
sults for the responses of a canard aircraft to step control deflec- 
tions and to the influence of two types of proportional automatic 
control. 

In all cases except the example illustrating the procedures of 
nonlinear mechanics ^ the nonlinear lariations have been ap- 
proximated by three linear segments. The principal means of 
studying the effiects of the nonlinearities were through the charac- 
teristics of the time responses in angle of attack and, in some 
cases, pitching velocity. The changes in the period and damping 
of the oscillations due to the nonlinearities have been discussed. 
The effiects of the autopilot proportionality constant on system 
stability were also investigated. The occurrence of continuous 
hunting oscillations was predicted and demonstrated for the at- 
titude stabilization system with proportional control for certain 
nonlinear pitching-moment variations and autopilot adjust- 
ments. 

INTRODUCTION 

In the classical treatments of the study of aircraft dynamics 
an assumption is usually made that the forces and moments 
resulting from small disturbances in position, velocity, and 
acceleration are linear variations with these quantities. 
However, some aircraft configurations do not exhibit linear 
force and moment characteristics in some flight conditions 
and, therefore, the linearized methods cannot be a satisfac- 
tory approximation over an appreciable range of the param- 
eter involved. Tliis report is concerned with the problem 
of aircraft dynamics as affected by nonlinear stability de- 
rivatives and, in particular, with the short-period longitudi- 
nal mode of motion of an aircraft flying at constant velocity. 
Consideration is also given to the problem of applying simple 
forms of automatic control to such an aircraft. 

The present work is presented in two independent parts. 
The first part considers several analytical approaches that 
may be used in the analysis and the study of the transient 


motions of the aircraft. These methods include the use of 
the operational calculus (Laplace transformation), the sta- 
bility criterion of Routh, and the theories of nonlinear 
mechanics. The second part presents some of the results 
obtained by using electrical analog equipment. The cases 
presented are for a hypothetical canard aircraft assumed to 
be flying at supersonic speed. 


SYMBOLS 
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”** ^ ac 
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dd 


Cr.= 


d6 


angle of attack, positive when nose is above 
relative wind vector, radians unless other- 
wise noted 


aircraft attitude or pitch angle, positive when 
nose is above horizontal, radians 


aircraft elevator deflection angle, positive 
when trailing edge is down, radians unless 
otherwise noted 
aircraft velocity, fps 
dynamic pressure, Ib/sq ft 
wing (reference) area, sq ft 
mean aerodynamic chord, ft 
mass of aircraft, slugs 
pitching moment of inertia, slug-ft^ 
differential operator, dfdt 
time 

pitching-moment coefficient. 

Pitching moment/giSc 
lift coefficient, Lift/gS 


> Supersedes recently declassified NACA RM L50L11, 1951. 
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h 

h 

a, 6, Cy d 

ay 7] 

f 

s 

K 


K, 

ai — a 


slope of a linear segment representing a por- 
tion of nonlinear Cm(«) function, per radian 
slope of a linear segment representing a por- 
tion of nonlinear Cl(o) function, per radian 
constants used in equations (6) and (9) and 
defined in appendix 

constants used in equations (2) and (3), 
respectively 

undamped natural frequency of aircraft alone, 
radians/sec 

nondimensional damping ratio of aircraft 
mode of motion 
Laplace transform variable 
general autopilot constant, ratio of control 
surface deflection to applied error angle 



Ui 


Si 


COy Cl 


X 


reference or desired angle of attack for closed - 
loop system with angle-of-attack feedback, 
radians 

reference or desired attitude angle for closed- 
loop attitude-stabilization system, radians 
constants required in cubic representation of 
nonlinear Cm(oc) function (eq. (12)) 
slope of trajectory in phase-plane example, 
dSlda 


ANALYSIS 

The analysis of the effects of nonlinear stability derivatives 
on the longitudinal motions of a canard aircraft is presented 
in three sections: the first is a statement of the problem and 
the assumptions involved; the second is a theoretical or 
analytical solution for the required transients and the study 
of stability; and the third is a discussion of the methods of 
nonlinear mechanics for these d}uiamical studies. 

STATEMENT OF PROBLEM 

The purpose of this report is to presen t^methods for deter- 
mining the nature of the transient responses of aircraft having 
nonlinear stability derivatives and to show some qualitative 
results of preliminary investigations. 

In the present report, consideration is given to the longitud- 
inal motions of a canard configuration with the angle-of- 
attack transients resulting from step elevator deflections 
serving to illustrate the effects of the nonlinear stability 
derivatives. Two longitudinal-stability derivatives were con- 
sidered to have nonlinear characteristics which were caused 
by the variations of pitching-moment coefficient and lift co- 
efficient with angle of attack. The main emphasis herein 
was placed on the effects of the nonlinear variation of pitch- 
ing-moment coefficient with angle of attack since the slope 
of this curve is a measure of the static stability of the aircraft. 
Figure 1 illustrates the nature of the nonlinearities under dis- 
cussion. Such nonlinearities would be apparent from wind- 
tunnel tests. 

In the analysis and study of the aircraft longitudinal 
motions, two degrees of freedom with constant forward 
velocity and disturbances from level flight have been as- 


sumed. This assumption results in a consideration of the 
short-period mode of motion only; therefore, only the effects 
of the nonlinearities on this mode will be apparent. 

In addition to studying transients of the aircraft in re- 
sponse to step elevator inputs, consideration was also given 
to the response characteristics resulting when such an air- 



Figure 1. — Typical nonlinearities in pitching-moment and lift 

coefficients. 




(a/ - a) 




Aircraft 


(a) 



(b) 


(a) Angle-of-attack stabilization, 

(b) Attitude stabilization. 

Figure 2. — Block diagrams illustrating the two types of proportional 
automatic control under consideration. 
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craft is subjected to simple forms of automatic stabilization. 
Two cases were investigated in this respect and are illustrated 
by block diagrams in figure 2. In both cases the automatic 
control consists of an error-sensing device and a zero-lag 
proportional servomotor. In figure 2 (a) angle-of-attack 
stabilization is considered with the control-surface deflec- 
tion 8 being proportional to the difference between a reference 
angle of attack at and the actual angle of attack a. In 
figure 2(b) attitude stabilization is considered where the 
control-surface deflection 8 is proportional to the difference 
between a reference attitude angle and the actual attitude 6. 


THEORETICAL METHODS 


Equations of motion. — The longitudinal equations of mo- 
tion of an aircraft having two degrees of freedom with 
constant forward velocity and disturbances from level flight 


in terms of the differential operator 


are 




mV 

~qS 


De- 


D9-C„{cx)-C^.^Da=C„,8 

^Da-C,ia)=C,,8 


( 1 ) 


The terms which are nonlinear functions of angle of attack 
are written as Cm{cc) and C^ia) and represent, respectively, 
the pitching-moment coefficient and lift coefficient that exist 
at the instantaneous angle of attack. 

Based on the usual assumptions of linear variations of Cm 
and Cl with a, equations (1) remain as ordinary linear differ- 
ential equations with constant coefficients and are handled 
by familiar mathematical methods. The types of non- 
linearities illustrated in figure 1, however, obviously make 
linearization of these parameters over the entire range of 
angle of attack an inadequate approximation. Since the 
methods for the solution of linear differential equations are 
so well known, an approximation of the nonlinear functions 
with linear segments, for which the equations can be written 
for a specified range of angle of attack, seems advisable. 
Figure 3 is an illustration of this method of representing a 
nonlinear function. The three straight lines shown are 
used to represent a variation of pitching-moment coefficient 
with angle of attack similar to that shown in figure 1 where 
the slope at outer a is twice the slope at central a. Each 
segment of this representation can be expressed as follows; 

Cm{cc)=kia+(T 


where ki and a must be specified for each range of angle of 
attack. In general, therefore, the representation of the 
nonlinearities are written as 


Gm(oc) = kia-\-a ( 2 ) 

Cl(«)=^2«+^ (3) 

where the values of kij k^j (7, and rj are constants and apply 
only for a specified range of angle of attack. 

The substitution of equations (2) and (3) into equations 
(1) results in linear differential equations with constant co- 
efficients, namely, 


^D«=C.,5+,' 


qS 


mV 


( 4 ) 


Da — k2a^CLi^Vv 


These equations apply only to the region of angle of attack 
in which the given constants are as defined. 

Aircraft transients. — Since variations in angle of attack 
following step elevator deflections may be in regions where 
the definition of the nonlinear functions has changed, com- 
plete knowledge of the motion is required and solutions of 
the equations must be written in terms of necessary condi- 
tions existing at an arbitrary time. In other words, the 
values of the angles and their rates of change must be known 
at the time one set of equations ceases to apply so that the 
same conditions and, therefore, continuity may be satisfied 
when the next set of equations becomes applicable. With 
these considerations the use of operational methods based 
on the Laplace transformation affords a means of handling 
the difficulties already noted and provides the required 
analytical expressions for the desired transient responses. 
References 1 to 4 give adequate discussions of the use of the 
Laplace transformation in the solution of differential equa- 
tions. Since the usual forms of the Laplace transformation 
require the initial conditions at zero time, it is convenient 
to make solutions for the required motions for each range 
of angle of attack and it may be necessary to redefine the 
time as zero whenever the equations must be changed be- 
cause of the representation of the nonlinearities. Such a 
procedure is obviously time consuming and tedious since 
solutions for all the parameters are required for a complete 
definition of initial conditions. A rather complete discussion 
of this procedure can be found in reference 5. 



Figure 3. — Illustration of the method of approximating nonlinear 
functions with linear segments. Cm(a) =A:ia4-a-. 
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Closed-loop stability. — The stability of a complete closed- 
loop system consisting of the aircraft and automatic-control 
unit, such as illustrated in figure 2, has been theorized by 
many authors when the system is linear throughout. (See, 
for example, refs. 2, 6, 7, and 8.) Briefly, this procedure is 
to plot the frequency response (ref. 9) of the aircraft as a 
Nyquist diagram (refs. 2, 7, and 8) or as an amplitude-phase 
plot (refs. 2 and 6) and then to determine the autopilot con- 
stant K required to make the system stable and give, as 
closely as possible, the desired response. Since the fre- 
quency-response method requires a completely linear system, 
such a method cannot be applied in the conventional manner 
and have meaning because of the definition of the nonlinear 
functions. 

This stability of the closed-loop system, however, can be 
examined by the methods of Routh and Hurwitz described 
in references 10, 11, and 12. This procedure requires the 
characteristic equations that exist for each region used in 
the representation of the nonlinear functions. If these in- 
vestigations indicate stability in all regions of the parametric 
representation, the stability of the closed loop is generally 
assured. Such methods could also be used, of course, to 
determine the values of the autopilot constant K required 
to assure complete stability. 

The application of this technique for the cases of angle-of- 
attack feedback and attitude stabilization is as follows: 

(a) Angle-of-attack feedback: The control equation is 
written as follows for the angle-of-attack feedback system 
(see fig. 2(a)): 

b=K\{ai — a) (5) 

When equation (5) is substituted into equations (4) and the 
Laplace transformations of the resulting expressions are 
formed, the solution for a can be written in terms of the 
transform variable as 

Ki{as^+bs)atis)+[l.C.(s)]a , 

where the term [I.C. (s)]a is the function of s (involving 
the initial conditions and constants) that results when the 
control equation and equations (4) have been transformed 
and solved for a(s) and the other coefficients are constants. 
These expressions are tabulated in terms of the aircraft 
parameters in the appendix. 

The stability criteria of the characteristic equation of 
equation (6) are 

2to}^+Kia>0 
(x)n-\-Kib^0 

The inequalities (7) allow a thorough investigation of the 
effects of Ki on stability. These criteria must be checked 
in each range of a used in defining the nonlinear functions 
which appear in equations (1) if the type of motion is to be 
completely investigated. The types of motion possible are 
similar to those of a simple spring-mass system with viscous 
damping, which also has a quadratic characteristic equation 
(see, for example, refs. 2 and 13). The equation which 
describes the attitude variations as this type of system 
operates can be determined by solving for 6(s) after equations 
(4) and (5) are transformed as described previously. 


The angle-of-attack feedback system does not exhibit a 
zero error at steady state even for the completely linear case. 
This condition may be seen by applying the final-value 
theorem (see, for example, ref. 2, pp. 72-73) to equation (6) 
and noting that the steady-state angle of attack does not 
equal the input at and, hence, the error at— a does not tend 
to zero. This characteristic also is present when the non- 
linear functions are considered. For the nonlinearities 
considered herein, all a transients are completely stable, that 
is, exhibit steady-state values. 

(b) Attitude stabilization: The control equation is written 
as follows when attitude stabilization is considered (see fig. 
2(b)): 

b=K2{Si-e) ( 8 ) 

When equation (8) is substituted into equations (4) and the 
transformations are made, the complete closed-loop solution 
for d{s) becomes 

. . K 2 (cs-]-d)di(s)-\-[I.C.(s)]e /nv 

S^ + 2tc.nS^ + (c^n^+K2C)s+K2d 

where [I.C.( 5 )]fl is the function of 5 (involving initial con- 
ditions and constants) that results when the control equation 
and equations (4) have been transformed and solved for 
6(s)y and the other factors are constants as defined in the 
appendix. 

The criteria for stability of the cubic characteristic equa- 
tion of equation (9) require that all coefficients be positive 
and that 

(2fco„) (o)„^+i^2^) — K2(i^0 (10) 

These criteria allow a complete investigation of the effect of 
K 2 on stability and should be checked for each angle-of- 
attack range where definition of the nonlinear functions 
changes. The angle-of-attack variation a(s) for this system 
can be solved for from the transformed equations (4) after 
the control equation (8) is considered. 

For the completely linear case when stability is assured, 
equation (9) reveals that the system is a zero-error system, 
the output tending to equal the input and the error (di—d) 
tending to equal zero. For the nonlinearities considered 
herein, two conditions affecting operation of the system can 
exist, and these conditions result primarily from the non- 
linear Cm{a) function. (For normal Ci,(q:) functions where 
the slope remains positive and is not greatly changed in the 
angle-of-attack regions, this nonlinearity will not affect the 
stability or oscillatory characteristics to any appreciable 
extent.) For the system having a nonlinear Cm(a) function 
and a K 2 satisfying the stability criteria for all values of a, it 
may be shown from the examination of steady-state a and 
steady-state 6 following a step input 6i that the system re- 
mains a zero-error system. Consider secondly the case 
where the Cm(oi) function has an unstable slope at central a 
and a stable slope at outer a and a K 2 such that stability is 
assured at outer a but not at central a. The instability at 
central a means that the angle-of-attack motion in this 
region is divergent to outer «; that is, no steady-state a can 
exist at central a. Upon examination of a{s) it will also be 
noted that no steady-state a can exist at outer a; however, 
the motion at outer a is stable, and the a tends to the region 
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of central « where the motion has been shown to be divergent. 
This process is summarized as follows: 

(1) At central a, the angle-of-attack motion is always 
divergent to outer a. 

(2) No steady-state a exists at outer a. 

(3) The tendency of the motion at outer a is toward the 
region of central a. 

It is concluded, therefore, that the angle-of-attack motion 
should, after some initial transients, become some form of 
continuous hunting oscillation (nonsinusoidal) rendering the 
complete system continuously oscillatory. A constant- 
amplitude, constant-period motion would also be suggested 
from the summary of the angle-of-attack tendencies when 
examined on an energy basis, an energy balance becoming 
established between that energy fed in by the unstable 
mode and that dissipated in the stable regions. 

In summary, the attitude-stabilization analysis reveals 
the following significant facts: 

(1) The autopilot proportionality constant K 2 can be 
adjusted for complete stability at all values of a and renders 
the system a zero-position-error system for the nonlinearities 
considered herein. 

(2) If the Cm{a) nonlinearity is such that the constant K 2 
assures stability at outer a but not at central a, continuous 
hunting oscillations can exist. 

Closed-loop transients. — Equations (6) and (9) are com- 
pletely expressed when ai{s) and Bi{s) are known. For a step 
function of unit magnitude, 

<9i(s) = ai(s)=7 

o 


this phase plane. The coordinates of the phase plane will 
usually allow the trajectories to describe successive states 
that the system may take as time proceeds but will usually 
not indicate the times at which the various states exist. For 
example, trajectories of displacement against rate of displace- 
ment would be sufficient to describe the successive states of 
motion of an oscillating system having one degree of freedom. 
The graphical representation of trajectories in the phase 
plane is facilitated by the construction of a family of isoclines 
where each isocline is a locus of points where the trajectories 
have a given slope (ref. 15, pp. 10-11 and 248-252). Inves- 
tigation of the singularities in the phase plane indicates the 
type of motions which result from given initial conditions, 
such as steady hunting oscillations or aperiodic or oscillatory 
responses that are either divergent or convergent. Refer- 
ences 16 and 17 are two examples of applications of the 
techniques of nonlinear mechanics to cases in which an 
automatic-control component has rendered a system essen- 
tially nonlinear. 

The example chosen for this discussion is an aircraft 
having a nonlinear variation of pitching-moment coefficient 
with angle of attack under the influence of angle-of-attack 
Stabilization (fig. 2(a)). The following assumptions are 
made: (1) the other coefficients are constants as listed in 
the appendix, (2) the reference angle of attack a, equals 
zero, and (3) the nonlinear pitching-moment variation has 
a positive slope for small positive and negative values of 
angle of attack. As a consequence of the second assump- 
tion, the following equation will apply: 

5= — K\a (11) 


Upon substitution into the proper equations, the transform 
function is known and the inverse transformation to a{t) 
or 0{t) or both ot{t) and Bit) can be determined. 

METHODS OF NONLINEAR MECHANICS 

For the types of nonlinear functions represented in figure 
1, the use of a cubic equation is suggested as a possible 
representation of the nonlinearity over an adequate range 
of angle of attack. Such a cubic representation would, of 
course, make the differential equations of motion of the 
aircraft nonlinear. To many such problems, however, the 
methods of nonlinear mechanics have proved a useful means 
of satisfactory analysis, although the results are seldom 
available as the familiar time responses. References 14 and 
15 afford very complete discussions of the methods and 
techniques used in applications of the theory of nonlinear 
mechanics. 

For the present report a single application of the methods 
of nonlinear mechanics will be discussed in detail to illustrate 
the application and usefulness of such methods to the problem 
considered. For the case chosen and, in general, for the 
cases where the principles of nonlinear mechanics are readily 
applicable, the procedure is to rewrite the differential equa- 
tions of motion into an independent system, that is, into a 
system of equations in which time appears only as a differ- 
ential dt. Such a representation is possible for many cases 
or can be derived by a suitable change in variable. After 
such a revision, a phase plane can be defined and the charac- 
teristics of the motion can be determined by trajectories in 


For the representation of the nonlinear function Cm(o)y 
consider the Cm curve in figure 1 having a positive slope at 
small value of a. If the a intercepts are symmetrical about 
a=0j a cubic representation of this curve would take the 
form 

Cm(o)=Coa^—Cia ( 12 ) 

where Co and Ci are constants. 

Considering a simplified notation for the equations of 
motion (eq. 4) and substituting equations (11) and (12) 
result in 

aiD^B-{-02DB—cooc^-{-Cia=a5(—Kia) 

biDB)- boDoc^r 630^=0 



where and Cm^ have been assumed to be zero. The com- 
plete definition of constants is given in the appendix. Equa- 
tions (13) are rewritten as follows to illustrate that they 
can be reduced to an independent system: 


dt a\ ai \ di cbi / 


(14a) 


n I 

*-*+" 57 “ 


dB 


where ^=-77 and use is made of the equation 62= ^61. 
at 


(14b) 

These 


equations suggest that the complete motions can be exam- 
ined by the instantaneous states of a and B. Tliis idea is 
also apparent from the Laplace transform procedure sug- 
gested in a previous section since only the initial values 
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of a and ^ are required for a definition of the motion. 
Therefore, the plane is the phase plane for these con- 
siderations. 

Before actual trajectories are drawn in this plane, a 
family of isoclines is determined. The slope of any path 
in the phase plane is dB\da and is determined by the division 
of equation (14a) by equation (14b). The resulting equa- 
tion is 


d^ 

da 


0/\ d\ 



ai 



I a 

di J 


(15) 


If dQjda is assigned a constant value X, equation (15) can be 
reduced to the form 


Co 





ai 



)a 

(h / 


(16) 


which is the isocline representing the locus of all points in 
the phase plane where trajectories must have the slope X. 

The singularities of the system are defined as the equilib- 
rium positions (steady-state values) that the system may 
attain and are determined from equations (14) by using 


M_da_ 

di~r^ 


(17) 


When the numerical values for the constants as given in 
the appendix are used, the following equation results for 
the isoclines: 

._-309(10)3a3+(4.5ix+257)a 

X+4.02 

where X is the slope in radians per second per radian. The 
singularities are found from equations (14) when conditions 
of equation (17) are imposed and are 

(1) «=0, 

(2) a=0.0278 radian= 1.59°, ^=0.1254 radians per second 

(3) a=— 0.0278 radian= — 1.59°, 0=— 0.1254 radians per 
second 

Figure 4 is a plot of the a,B phase plane and shows some 
isoclines for various X. It is important to note that the 
trajectories or paths which describe the motion have a 
definite direction that depends on the position in the phase 
plane. This direction is determined from equations (14) 
by considering the value of equation (14a) to be a vector 
component in the d direction and the value of equation (14b) 
to be a vector component in the a direction. The direction 
of the resultant of these two components is the direction of 
the trajectory at that point in the phase plane; this resultant 
vector is termed the phase velocity at that point since the 
vector components determined from equations (14) are them- 
selves velocities. The arrows shown on the isoclines are not 
phase-velocity vectors but indicate only the direction of the 
phase velocity on the various isoclines. 



Figure 4. — The phase plane (a, e) for the example illustrating the methods of nonlinear mechanics. Three trajectories are shown; isoclines of 

constant X are labeled. 
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Several trajectories are shown in figure 4 for various initial 
conditions; these were drawn by the graphical method dis- 
cussed in reference 15, pages 248-252. It is evident from 
these trajectories and from the directions of phase velocity 
on the given isoclines that the two singularities (dzl.59°, 
±0.1254) are points of stable equilibrium; that is, the 
motion following a disturbance will oscillate about and tend 
to one of these stable states. A singularity of this type is 
termed (refs. 15 and 16) a ^‘stable focus"' which means that 
the motion tends to this equilibrium in an oscillatory manner. 
The other singularity at (0,0) is an unstable equilibrium 
point. This instability is evidenced by the fact that the 
slightest disturbance from this point results in a motion to 
one of the stable foci. This singularity at the origin is 
termed a ‘^saddle point." 

In tills brief discussion it is evident that much can be 
learned about the characteristics of the system without 
resorting to actual solution of the nonlinear differential 
equations. The results presented for this case from an 
inspection of the nature of the phase plane can also be 
verified by the application of the various theorems presented 
in references 14 and 15, which allow a determination of the 
stability of singularities and of such an important property 
as the existence or nonexistence of steady oscillations in the 
system. The methods of nonlinear mechanics, however, are 
relatively new and, in their present state of development, 
they cannot be expected to handle an extreme variety of 
problems. The methods become quite involved even for 
systems having two degrees of freedom in which definition of 
the motion requires more than two quantities; the methods 
become almost impossible when more degrees of freedom are 
introduced. 

ANALOG RESULTS AND DISCUSSION 

As a part of the investigation of the effects of nonlinear 
stability derivatives on the transient longitudinal motions 
of an aircraft, the problem was studied through the use of 
electrical analog equipment. The analogs employed were 
the Engine and Control Simulator manufactured by Phil- 
brick Researches, Inc., and the Reeves Electronic Analog 
Computer (hereinafter designated the Philbrick Analog and 
the RE AC, respectively), the units being, respectively, at the 
Lewis Flight Propulsion Laboratory and the Ames Aero- 
nautical Laboratory of the NACA. Descriptions of the 
analog units are available from various sources, for example, 
references 18 and 19. 

The nonlinearities considered in these studies included 
both nonlinear functions represented in equations (4), that is, 
the nonlinear variations of pitching-moment coefficient and 
lift coefficient with angle of attack. In all problems pre- 
sented for both analogs, the nonlinear functions were repre- 
sented by straight-line segments although the methods of 
representation used in the two analogs differed. The other 
basic constants used are tabulated in the appendix, and the 
nonlinear functions arc defined when presented. The 
presentation and discussion of results are given in two sec- 
tions: (1) aircraft transient responses and (2) closed-loop 
transient responses. The first section contains analog solu- 
tions of the transients in angle of attack (and, in some cases, 
rate of pitching) resulting from step deflections of the canard 
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elevator control surfaces. The second section presents some 
solutions for the closed-loop cases shown in figure 2. In both 
sections various types of nonlinear functions were chosen as 
being possibly representative of the nonlinearities likely for 
the assumed canard configuration. 

The results from the Philbrick Analog were obtained by 
photographing an oscilloscope; the REAC results were made 
with recording elements. The intent of the Philbrick studies 
was to survey qualitatively a wide range of nonlinear func- 
tions by observing and recording the more unusual tran- 
sient conditions. For this reason high accuracy in aircraft 
simulation was not required and the results are discussed in 
this light. In the REAC problems more exact solutions were 
required, the accuracy in final a and d transients being esti- 
mated at approximately 2 percent and the accuracy of 
individual REAC components being considerably better. 

AIRCRAFT TRANSIENT RESPONSES 

The analog results presented herein are aircraft transient 
responses in angle of attack (and pitching velocity) due to 
step deflections of the canard elevator surfaces. Although 
only a few nonlinear situations are discussed in detail, the 
remarks which summarize the report are based on a much 
larger number of cases for which analog solutions were ob- 
tained. 

Case I. Nonlinear Cm{a): 0<n< 1, — 2°<a<2°; 

dC 

a>2®. — In case I the nonlinear function 

is the pitching-moment coefficient due to angle of 
attack where the slope at values of a less than —2° 
and greater than 2° is larger than the value of the slope in the 
region — 2°<a<2°. In this case, the slope is always nega- 
tive and stability is assured. Since the frequency of the 
transient oscillations is determined principally by this slope, 
changes in the frequency (period) are e.xpected as the region 
in which the value of the steady-state angle of attack changes 
because of various elevator deflections. Little change in the 
damping envelope is expected. 

Results from the Philbrick Analog are presented in figure 
5. In this run the value of ki was equal to —5.14 and the 
slope for — 2°<a<2° is 0.463A:i= — 2.38. The a transient 
shown in figure 5(a) is the completely linear case where 
Ctn(oc)=kia. The total transient time shown represents 
1.325 seconds. When the described nonlinearity in Cm(oc) 
is introduced, the results are shown by the four remaining 
records of figure 5. Curves of figures 5(b) and 5(d) are 
two a transients, whereas curves of figures 5(c) and 5(e) are 
the corresponding nonlinearities. These nonlinearities were 
photographed by using the output of the element generating 
the nonlinear function as the vertical coordinate of the 
oscilloscope and the angle of attack as the horizontal co- 
ordinate. Hence, the region of angle of attack covered in 
the transients of figures 5(b) and 5(d) is shown again by the 
region of the nonlinearity pictured. The large spots shown 
on the nonlinearities represent the trim values of a due to the 
old and new 8 input. Hence, figure 5(c) shows that the 
values of 8 change in such a way as to make the trim values 
of a change from one extreme of the nonlinearity to the other. 
Since the slopes at these outer values of a are the same as that 
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(a) Linear case, -^= — 5.14. 

(b) a transient. (c) Nonlinearity, 

(d) a transient. (e) Nonlinearity. 


Figure 5. — Philbrick solutions for case I. 

f \ 


dCm 


Aircraft a transients and parts of nonlinearity involved; nonlinear Cm(«), 


dC, 


dCL 


> 20 . !1^=-2.38, -2°<a<2°; ^^=3.49, all a. 


da 


shown for the linear case, the period should remain prac- 
tically unchanged (compare figs. 5(a) and 5(b)). 

Figure 5 (d) is the a transient when the change in 8 causes 
the trim a to ^change from a point where the slope is large to 
one in the region of reduced slope. The ratio increase in 
period (reduction in frequency) shown in figure 5 (d) is ap- 
proximately equal to the theoretical variation of being in- 
versely proportional to the square root of the decrease in 
slope. (This variation will apply for lightly damped sys- 
tems.) For the three transients shown, the damping envelope 
is seen to be essentially unchanged, the time required to reach 
steady state being about 1 second in each case. 


Two solutions from the REAC for this case (reproductions 
of actual recordings) are presented in figure 6. For these 
solutions the value of ki was —6.0 and n=0.5, and the cases 
shown are the results of moving the elevator in a programmed 
manner from 0° to about 8° and after steady state is reached 
to about —4°. The variations of 6, a, and 8 with time are 
shown in tliis figure. The slight phase difference between 
control motion and the transients in tli6 REAC solutions are 
due to the recording equipment and not to aircraft dynamics. 
The ±2° (dzO.035 radian) positions noting the break points 
in the nonlinearity are also indicated on the a transient. The 
same general characteristics are observed for these two tran- 
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Figure 6. — REAC solutions for case I. Aircraft transients; nonlinear Cm(a) function with ^^= — 3.0, 

-2°<a<2°; ^=3.49, all a. 


sients as mentioned previously, namely, the change in period 
and approximately the same time required for both transients 
to reach steady state. The steady-state values for a. and d 
check reasonably well with those predicted theoretically. In 
figure 6 the responses show clearly that when the angle-of- 
attack variations are such as to be in the region of breaks in 
the nonlinearity, the periodicity is not constant; when they, 
however, are in an a region where the slope of the nonlinearity 
is constant (near steady state), the period does remain con- 
stant. For cases such as this illustration where it is rela- 
tively certain that the system is very lightly damped, this 
change in periodicity and the eventual constant period near 
steady state can be attributed to a nonlinearity in the effec- 
tive spring constant of the system, in this case (a). 

Case n. Nonlinear «>2°; 

—nkiy 0<^7i<;i, — 2°<^a<^2°. — In case II the variation of the 
pitching-moment coefficient with angle of attack is nonlinear 
and considered to have a positive (unstable) slope in the 
range — 2°<Co:<.2°. A trim angle of attack for a given ele- 
vator deflection 6 cannot be expected to exist in the region 
— 2°< 2° since in this angle-of-attack region the transients 

are divergent. (This condition is shown by the presence of 
two real roots of opposite signs of the quadratic characteristic 
equation of equations (4), the negative root indicating a 
damped subsidence and the positive root indicating the diver- 
gent aperiodic mode.) 

Figure 7 presents some Philbrick solutions for this case 
where Z:i=— 0.838, n — l. The total transient time shown 
represents 2.13 seconds. The a transient shown in figure 7 (a) 


is for the completely" linear case with Cm{(x)=kia. The period 
of the damped oscillation is 0.284 second. Figures 7 (b), 
7 (d), and 7 (f) are a transients for different b conditions, and 
figures 7 (c), 7 (e), and 7 (g) are the parts of the nonlinear 
function involved in each transient. 

From the nonlinearity shown in figure 7 (c) the change in 
trim angle of attack is large with the new trim angle of at- 
tack well into the region of stable slope. Therefore, little 
change is expected or noted in the transient of figure 7 (b) 
as compared with the linear case in figure 7 (a). As the 
value of 5 is changed and the new trim angle of attack comes 
nearer to the break in the nonlinearity (see fig. 7 (e)), the 
initial portion of the transient of figure 7 (d) has become 
distorted. The characteristics of the oscillations about the 
new trim angle of attack, however, remain unchanged. As 
the b is changed only slightly, a new trim angle of attack is 
found to the left of the region of unstable slope as shown in 
figure 7 (g). The first overshoot of this transient (fig. 7 (f)) 
is to an angle of attack in the unstable region, but the 
energy of the system is such that the complete unstable 
region is not traversed and the oscillatory characteristics 
about the new trim position become established. 

Only a very slight change in the b was found to be required 
for the change in trim angle of attack exhibited in the 
transients of figures 7 (d) and 7 (f). Only when the new 
trim angle of attack was near the breaks in the nonlinearities 
did the transients indicate any peculiarities which might 
give an indication of the presence of a nonlinearity in the 
system. 
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(b) a transient. (c) Nonlinearity. 

(d) a transient. (e) Nonlinearity. 

. ^ (0 « transient. (g) Nonlinearity. ^ 

dCm 

Figure 7. — Philbrick solutions for case II. Aircraft a transients and parts of nonlinearity involved; nonlinear Cm (a), -^^= — 0.838, a< — 2°, 


«> 2°; -^=0.838, -2°<a<2” 
da 


; ^^=3.49, all a. 


The REAC solutions presented in figure 8 for this case 
are for a value of A:i=— 3.0, n=0.5. The two examples 
shown start from the same values of 5, trim a, and d with 
the first changing from to 8 2° and the second 

changing to 4°. In both of these conditions, the 

angle-of-attack variations are from a point of stable equi- 
librium at positive a to a stable trim position at negative a. 
The transients indicate that the initial portions give the 
only hint of nonlinearity in the system. As noted previ- 
ously, this initial part of the a transients becomes more out 


of the ordinary as the new position of trim a approaches the 
break point of the nonlinearity. These slight slope changes 
in the a transient are almost unnoticeable in the second 
curve. Examination of the 6 transients, however, indicates 
clearly the effect of the nonlinearity 'on the rate of pitching 
in the initial phase of the response. After this first effect 
on the transients, the oscillatory characteristics become very 
regular. 

rjn 

Case HI. Nonlinear Cl(o): 2°, a>2°; 
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dCL 

da 


= 3.49, all a. 


-^=n)t 2 , 0<w< 1,— 2°<a<2°.— In case III the variation of 
da 

lift coefRcient with angle of attack is the nonlinear function 
having a larger slope at the outer values of a. Since the 
aircraft is lightly damped, the primary contribution of the 
lift-curve slope to the oscillatory characteristics is in the 
damping of the motion. Since, for this case, the slope of 
the pitching-moment-coefRcient curve against angle of 

da 

in the period of the transient oscillations. 

In figure 9 three Philbrick solutions are shown. Figure 
9 (a) is the completely linear case with ko=3A9 (n=1.0) 
and A:i= — 1.90. The transient time is 2.48 seconds, the 
period is 0.194 second, and the cycles to damp to one-tenth 
amplitude are 3.39. Figure 9 (b) is the transient which 
results when the trim a changes from a region of large slope 
(dCJda) to a region of a of reduced slope. The value of n 
is 0.328 in this instance, that is, the slope at central a is 
about one-third of the slope at outer a. The period is 
essentially unchanged although the number of cycles to 
damp to one-tenth amplitude has increased to 4.20. For a 
change in trim a from the region of reduced slope to the 
region where the slope is the same as the completely linear 
case, the transient in figure 9 (c) is essentially the same as 
the one shown in figure 9 (a) for the linear case. 

The change in slope used in this example is much greater 
than would be expected in a practical case; thus this change 
indicates that the actual effects of this parameter might result 
in very little change in the a transient characteristics. It is 


also emphasized that a simple examination of the transient 
(or even a comparison between figs. 9 (b) and 9 (c)) could 
not indicate that the nonlinearity was in lift-curve slope 
because the effect of variation of damping in pitch 
would also be manifested by changing the damping but not 
affecting the period. In actual flight tests, however, normal- 
acceleration records might help to isolate this effect as a 
nonlinear Cl{ol). 

The RE AC solutions shown in figure 10 have A:2=4.19 and 
7i,=0.834. The linear equation of = — 3.0a was also 

used in these runs. The variations of 8 shown were from 
5=0° to 5«8° to 4°. The periodicity is constant for 

these cases even though the variations of a in the second run 
are about the break point (a= —0.035 radian) of the non- 
linearity. Since the two slopes of the nonlinear Cx(a) func- 
tion are nearly the same, the damping is practically unaffected 
in the two cases; in fact, plots of the logarithm of the peak 
amplitudes about steady-state values against time for these 
two cases indicate essentially the same slope and, therefore, 
the same damping, although the second run has a slight 
amount of scatter in the points. These cases, which have a 
more reasonable value for n than was used in the Pliilbrick 
solutions, indicate the slight effects of this parameter under 
the conditions investigated. 

Case IV. Nonlinear C^ia) and (7i,(a).— In case IV the two 
nonlinearities used in the REAC solutions for cases I and 
III were used together in the same run for two transient solu- 
tions. These results are shown in figure 11 for 5 changes of 
0° to 4° and 0° to 8°. In the first transient only the first 
overshoot is greater than a=2°, the break point of the non- 


attack is a constant 


( 


no change is to be expected 
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(a) n=1.0, completely linear case. 

(b) w= 0.328. 

(c) n= 0.328. 


Figure 9. — ^^Philbrick solutions for case III. Aircraft a transients; 
nonlinear Cl(«); ^=3.49, a<-2°, a>2°; ^=n(3.49), 

-2°<«<2°; ^=-1-90. all “• 

linearities; hence, the transient characteristics are those 
defined entirely by the slopes of the nonlinearities at small a. 
In the second transient the period becomes constant after 
the first two cycles and exhibits characteristics similar to 
that of the first transient of figure 6 with the single nonlinear 
function Cm{ot). These results might have been expected 
since the nonlinear Cz,(a) had so little effect in case III. 
Other runs for cases where two nonlinearities were in the 
system indicated that the most noticeable changes in the 
transients occurred with changes in the Cmiot) function, as 
indicated previously. 

Review of effects of nonlinearities on aircraft transients. — 
Nonlinear variations of the pitching moment resulted almost 


solely in effects on tlie period of the transient oscillations. 
Differences in the period of the oscillations in the different 
regions of angle of attack, due to the different slopes of the 
nonlinear function in these regions, were an effective way.^ 
noting the presence of this nonlinearity. Another noticeable 
effect was the irregularity in the period of the oscillations 
when the trim angle of attack for the various elevator deflec- 
tions was near a break point of the nonlinearity. Wlien the 
trim angle of attack is removed from the break in the non- 
linearity, the oscillation characteristics become very regular, 
that is, constant period and damping. If only moderate 
changes in trim angle of attack result from tl^ie control motion 
and the nonlinearity or a part thereof is traversed, the initial 
sections of the transients may indicate the presence of the 
nonlinearity through small slope variations, although the 
final oscillatory characteristics are quite regular. For large 
changes in trim angles of attack under the conditions noted, 
the small variations in the initial sections of the transient 
may become so slight as to be indistinguishable from the 
perfectly linear case. 

For cases when the nonlinear variation was in lift coeflS- 
cient, large slope changes were required before changes in the 
oscillatory characteristics became appreciable for trim angles 
of attack in the different regions. This nonlinearity affected 
only the damping characteristics of the motions, the period 
remaining constant. This effect could be expected since 
very large changes in damping would be required to alter the 
period of this lightly damped system. 

Wien the two nonlinearities were considered simultane- 
ously, with reasonable variations in each, the predominant 
changes in transient characteristics were similar to those 
noted for the single nonlinear pitching-moment variation, 
that is, the changes in periodicity previously mentioned with 
the time to damp remaining essentially unaltered. 

CLOSED-LOOP TRANSIENT RESPONSES 

The results in this section are transient response curves 
illustrating the dynamics of the two closed-loop systems 
depicted in figure 2. The discussion is made in two parts: 
the closed loop with angle-of-attack feedback and the atti- 
tude stabilization; Philbrick and REAC solutions for several 
cases are included. The nonlinear characteristics are those 
previously described for the aircraft; the error-sensing devices 
and proportional servomotors are considered as zero-lag 
systems. The transients are the responses to step input, 
command signals in at and 0,. 

Angle-of-attack feedback. — The block diagram for the sys- 
tem with angle-of-attack feedback is shown in figure 2 (a). 
Transients from the Philbrick Analog are ^hown for two 
cases: one where the nonlinear function is Cl{o) and the 
other for the nonlinearity Cm{ot)> REAC solutions were 
made for two different nonlinear funqtions of pitching- 
moment coefficient. 

In the Philbrick solutions the value of Ki was approx- 
imately unity. For the case with nonlinear Cl{ol)^ the 
conditions are the same as described for case III of the 
previous section where the lift-curve slope at outer values of 
a is greater than the slope at small values of a. The tran- 
sients shown in figure ^12 resulted from changes in oli that 
caused the final values of a to be in regions where the lift- 
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FiauRE 10. — HEAC solutions for case III. 


Aircraft transients: nonlinear C/,(a) function with ^^= 4 . 19, o;<[ — 2°, a>2°; -^^=3.49, —2°<^a<^2°; 

^’!=_3.0, all a. 
da 


curve slope is different. Figure 12 (a) shows the a transient 
about a steady-state value of a in the region where the lift- 

curve slope is reduced ^-^-^=1.14^* The oscillations in 

figure 12 (b) are about a steady-state value of a in an a 

region of larger lift-curve slope The increase 

in lift-curve slope is seen to improve the damping character- 
istics of the closed-loop system with angle-of-attack feed- 
back in much the same way as the aircraft transients were 
improved. Investigation of the characteristic equation of 
tliis closed loop, the denominator of equation (6), shows that 
the principal effect of such a change in lift-curve slope would 
be most evident in the damping. Comparison of transients 
of figure 12 as well as an examination of the characteristic 
equation indicated that the frequency of the oscillations is 
constant. 

The Philbrick solutions with nonlinear Cm(a) were made 
for the case where the slope of the pitching-moment-coeffici- 
ent curve was positive (unstable) at small a and stable at 
outer a. The outer slopes were —1.90, and the center slope 
was 0.655. Two a transients illustrating the effects of this 
nonlinearity on the closed-loop responses are given in figure 
13. In figure 13 (a) the change in reference angle at results 
in a change in the steady-state values of a from one region 
of stable slope to the other. In figure 13 (b) the input 
change causes the final value of a to be a value in the region 
of angle of attack where the aircraft has an unstable pitching- 


moment curve, and the transient is satisfactory. The 
physics of this change is explained by the action of the angle- 
of-attack feedback which causes the control-surface deflection 
8 to give a moment proportional to the instantaneous a; 
thus, this moment adds to the effective spring constant of 
the system. In this instance, the value of Ki is large enough 
so that tliis contribution overcomes the destabilizing effect 
of the spring in the aircraft at central a and renders the 
overall system stable in this region of a. This stability could 
also be determined from the inequalities (7) wliich show that 
the value of Ki used assures complete stability of this closed- 
loop system. 

REAC solutions were also made for a case where the pitch- 
ing-moment curve had an unstable slope at small values of 
angle of attack. The same conditions as discussed for 
REAC runs in case II were used in this solution with a value 

of ^^=—3.0 for the stable slopes at outer a and 

for — 2°<Ca<2°. The value of Ki was 1.0. The three a 
transients of figure 14 are for the at step inputs of 1°, 4°, and 
8°. The 2° position noting the break in the nonlinearity is 
shown on each curve. For this case, an examination of the 
inequalities (7) reveals that the value of Ki = l.O is not large 
enough to assure complete stability in the region — 2°<a<2^; 
hence, no steady-state value of a is possible in this region. 
Therefore, for even the small input 1°, the final value of a 
is greater than the break point of the nonlinearity. The 
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Figure 11. — REAC solutions for case IV. Aircraft transients; nonlinear Cm(a), ^^= — 6.0, a< — 2°, a>2°:^^ 

da da 

Cl(oc), ^=4.19, a<-2°, a>2°; ^=3.49, -2°<«<2°. 


= — 3.0, —2°<a<2°; nonlinear 



measured values of steady-state a due to the inputs involved 
are 2.49°, 3.21°, and 4.24° compared with the theoretical 
values of 2.50°, 3.26°, and 4.30°, respectively. Since the 
value of a is greater than 2° after the first cycle for all the 
transients shown, the period of the oscillations is essentially 
the same for all cases. 

The final REAC solutions are shown for the case where the 
Cm(ci) nonlinearity is the same as case I, namely, ^^= — 3.0 

fin 

where — 2°<a<^2° and when a<— 2°, a>2°. 

The inputs were again 1°, 4°, and 8°, and the value of 
was 1.0. The resulting closed-loop a transients are given in 
figure 15. In these runs too, the steady-state values of a 
check satisfactorily with the theoretical values. Note that 
figure 15 (c) for ai=8° has a final value for a of about 2°; 
hence, the oscillations are about the break point of the 
nonlinearity. The period of this transient is very constant 
and is diflferent from the period of figures 15 (a) and 15 (b) 
which both remain in the region 0<a<2°. The period of 
figures 15 (a) and 15 (b) is about 1.17 times the period of 
figure 15 (c). A theoretical value of the change in period, 
obtained from consideration of the values defined by the two 
slopes of is 1.15; this value is close to the measured 

value even though the oscillations were not in the region 
where this result could be expected to apply. Stability of 
the combination of the aircraft and autopilot with control 


proportional to angle of attack was obtained by adjustment 
of the autopilot proportionality constant, and the effects of 
the nonlinearities were similar to those noted for the aircraft 
transients to step control deflections. 

Examination of figures 14 and 15 indicates that all the 
cases shown with nonlinearities in the system exhibit a 
steady-state error between the reference and the actual a. 
This characteristic is not due primarily to the nonlinearities 
of the system. An investigation of a linear aircraft in this 
system of angle-of-attack feedback reveals that the complete 
system is not a zero-position-error system; therefore, unless 
some form of compensation is included in the loop the output 
will not equal the input at steady state. For the cases shown, 
however, adequate compensation of the system with non- 
linearities as described may be no more difficult than that 
required to make the linear system a zero-position-error 
system. 

Attitude stabilization. — In the results (REAC) to be pre- 
sented for the problem of attitude stabilization, the only 
nonlinearity considered was Cm{a). In the first example the 
Cm{oi) function used for REAC solutions for case I, dCmIda 
being negative for all a but having a larger slope at outer a, 
was employed. The second example is the Cm{oi) function 
for case II, dCJda having a stable (negative) slope at outer 
a and an unstable slope at — 2°<a<2° and a K 2 which does 
not assure stability at central a. The analysis in previous 
sections indicated that the first example should exhibit 
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dC 

(a) Steady-state a in region where -^=1.14. 


(b) Steady-state a in region where ^^=3.49. 

Figure 12. — Philbrick solutions for proportional control system with 

dC 

angle-of-attack feedback, a transients, nonlinear C^Ca), -^=3.40, 
a<-2°, a>2°;^=1.14, -2“<a<2°; - 1.90, all a; auto- 

pilot constant Ki « 1.0. 

constant steady-state values after transients die out, whereas 
the second example should show continuous hunting oscilla- 
tions. For both cases K 2 = 1 .0. 

Figure 16 illustrates the first example where the Cm{ot) has 
a negative slope at all a. Figure 16 (a) is for 0^=4. 6° and 
figure 16 (b) is for 6i=9^. The upper cuives in figure 16 (a) 
and figure 16 (b) are the a variations, whereas the lower 
curves are of the error df—d. In both of these cases the 
error tends to zero indicating the zero-position-error charac- 
teristic of the system. The cubic characteristic equation of 
this system, see equation (9), would show an oscillatory 
mode of motion superimposed upon the convergent aperiodic 
mode and is indicated by the negative real root. From the 
characteristics of the system at central a, the oscillatory 
mode damps to one-half amplitude in about one-third of 
the time required for the aperiodic mode to damp the same 
amount. In figure 16 (b), the period becomes constant only 
after the variations of a become less than 2®, whereas in 
figure 16 (a) the period is constant throughout. 

The results for the second example where the Cm(oi) func- 
tion had an unstable slope in the region — 2°<a<2° are 
presented in figure 17. Curves in figures 17 (a), 17 (b), and 


17 (c) were results of step inputs of df of 0.7°, 4.6°, and 8.7°, 
respectively. In each set the time variations of a and 6i—6 
are shown. Though the initial portions of the responses are 
different for the three cases, the resulting hunting oscillations 
are essentially the same. These curves show the constant- 
amplitude constant-period oscillations predicted in the 
theoretical analysis for this type of nonlinearity. 

Stability of the attitude-stabilization system, when the 
control-surface deflection is proportional to the error in 
attitude, could be completely assured for the types of non- 
linearities considered herein with proper adjustment of tlie 
autopilot proportionality constant. The existence of con- 
tinuous hunting oscillations for a particular nonlinear 
pitching-moment variation and improper autopilot adjust- 
ment was predicted theoretically and demonstrated by 
analog solutions. The pilot of a conventional aircraft might 
apply control in a manner that is approximately proportional, 
though probably with a time lag, and such a hunting oscilla- 
tion as described might be realized if that aircraft had a 
characteristic like the nonlinearity in the data of figure 17. 




(b) Steady-state a in region where ^^^=0.655. 

Figure 13. — Philbrick solutions for proportional control system with 
angle-of-attack feedback, a transients and parts of the nonlin- 

earity; nonlinear C„(a), «>2°; ^ = 0.655, 

J 

— 2°<a<2°; -^=3.49, all a; autopilot constant K\ « 1.0. 
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a, radian 


0 - 


. 10 ^ 



(b) 



(a) ai=l\ 

(b) = 

(c) «i=8°. 

dC 

Figure 14. — REAC solutions for proportional control system with angle-of-attack feedback. Nonlinear C„,(a), — 2°<a<2°; 

^=-3.0, a<-2°, a>2®; ^=3.49, all a; autopilot constant ^ 1 = 1.0. 
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(a) 


a, radian 


0 


.04 




(b) a, =4°. 

(c) 

Figure 15. — REAC solutions for proportional control system with angle-of-attack feedback. Nonlinear C,„(a), — 2°<q:<2°; 


dCm 

da 


6.0, a<-2°,a>2°; 


dCL 

da 


3.49, all a \ autopilot constant 




1 . 0 , 
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(a) 0i = 4.6°. 

(b) 0,=9'’. 


Figure 16. — REAC solutions for proportional control for attitude stabilization. Nonlinear C„(a), 


a>2°; ^^^=3.49; autopilot constant fv 2 = 1 . 0 . 


dC„ 

da 


3.0,-2°<a<2°; 


dC„ 

da 


-6.0, a<-2°. 
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(c) Time, sec 


(a) e 

(b) e 

(c) e 


=0.7°. 
= 4.6°. 
= 8.7°. 


Fioure 17. — REAC solutions for proportional control for attitude stabilization. Nonlinear CmM, ^^=1.5, — 2°<'a<^2°; ^^-=- 

da da 

dCu 


■3.0.a< -2° 


«>2°; -^"=3.49; autopilot constant K 2 —I.O, 
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CONCLUSIONS 

A study has been made concerning the effects of two non- 
linear stability derivatives, caused by the nonlinear varia- 
tions of pitching-moment and lift coefficients with angle of 
attack, on the longitudinal motions of an airciaft. In 
addition to the theoretical methods discussed, the excellent 
adaptation of analog equipment to this dynamical study 
gave results for the responses of a canard aircraft to step 
control deflections and to the influence of two types of 
proportional automatic control. From the theoretical con- 
siderations and the analog results presented herein, the 
following conclusions summarize the investigation: 

1. When the nonlinear functions are approximated by a 
series of linear segments, the methods of the Laplace trans- 
formation allow the transient responses of the aircraft alone 
and also with proportional automatic control to be calcu- 
lated. The computational procedure required, however, is 
tedious and time consuming. 

2. The most noticeable effects of the nonlinear pitching- 
moment variations on the response of the aircraft to step 
control deflections was in the periodicity of the oscillations. 
The occurrence of an irregular period and of different 
constant-period oscillations about the various trim angles 
of attack was the most satisfactory way of noting the presence 
of this nonlinearity. 

3. Nonlinearity in the lift curve affected only the damping 
of the transient angle-of-attack oscillations, the period 
remaining constant. Large slope changes in this nonlinear 
function were required before the effect on the time to damp 
became appreciable. 

4. Simultaneous occurrence of the two types of nonlinear- 
ites considered herein, with reasonable values of each, resulted 
in transient characteristics which were similar to those with 
the single nonlinear pitching-moment function. 

5. Stability of the combination of the aircraft and auto- 
pilot with control proportional to angle of attack was ob- 
tained by adjustment of the autopilot proportionality con- 
stant, and the effects of the nonlinearities were similar to those 
noted for the aircraft transients to step control deflections. 

6. Stability of the attitude-stabilization system, when the 
control-surface deflection is proportional to the error in 
attitude, could be completely assured for the types of non- 
linearities considered herein with proper adjustment of the 
autopilot proportionality constant. The existence of con- 
tinuous hunting oscillations for a particular nonlinear 
pitching-moment variation and improper autopilot adjust- 
ment was predicted theoretically and demonstrated by analog 
solutions. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., December 29^ 1950. 


APPENDIX 


FACTORS AND COEFFICIENTS USED IN THE ANALYSIS 

The complete expressions for the factors and coefficients 
used in the analysis section of the report, in the order in 
which they appear in the text are as follows: 



^_ a2&4 I 0-5 

cb\h2 eb\ 


[I.C. (s)]„=a(0)s+[^^‘ a(0)+^(0)+^J 


d\ Cb\h2 


cyy. 02,(l^,tC2 
2rC0n=— -1-— 

dl d\ 0i 


2 0/JC2 

COn = 7 

di02 dl 

(When and 2fco;^>0, corresponds to the undamped 

natural frequency of the short-period mode of motion of the 
aircraft alone and f represents the conventional nondimen- 
sional damping ratio.) 

d^ (X464 
dl dib2 


ajco I kib^, 

dib2 d\b2 


[I.C. (s)],=0 (0) (O)(|+^+0+e (0)]s 


+ 




+ 


b2(T—d^r] ^ 1 


dib2 


di= 


dib2 

Iz. 

qSc 


) 




C 


Cti — Cn- 


2V 


ai=C, 


mf 






Li 
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Values used for example of method of nonlinear mechanics 
are as follows: 

Co=-546 
Ci= — 1.5 
ai=0.00177 
ao=0.00712 

Cl4 — 0 

as= 1.045 

-62=0.774 

64=0 



Xi=1.0 

For the solutions made on the Philbrick and Reeves 
Analog equipment, the numerical values for the afore- 
mentioned constants Ui, a 2 , a4, 61, 62 , and 64 were used. 

The values for these coefficients resulted from the following 
flight conditions and stability derivatives: 


Mach number= 1.8 
F=2009 ft/sec 
^=4800 Ib/sq ft 
S=2.52 sq ft 
c = 1.4 ft 


oc ~ 
2V 


-20.43 


n 045 

b8 ^ ^ 




/y=30 slug-ft^ 
m=4.66 slugs 
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